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We measure the transverse thermoelectric coefBcient a^y in simulations of type-II superconductors 
in the vortex liquid regime, using the time-dependent Ginzburg-Landau (TDGL) equation with 
thermal noise. Our results are in reasonably good quantitative agreement with experimental data 
on cuprate samples, suggesting that this simple model contains much of the physics behind the large 
Nernst effect due to superconducting fluctuations observed in these materials. 
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I. INTRODUCTION 

The Nernst effect in the cuprate superconductors has 
recently become a focus of attention both experimen- 
tally Q, S H H H and theoretically H 0, B 11. The 
Nernst effect is the electric field induced when the sam- 
ple has a temperature gradient, VT, perpendicular to 
the magnetic field, H; this electric field is perpendicular 
to both VT and H. For some cuprate superconductors 
the Nernst effect due to superconducting fluctuations is 
detectable at ternperatures far above the transition tem- 
perature, 0, 0. Ussishkin et al. ,7] calculated the 
low-field Nernst effect for T > due to superconducting 
fluctuations, obtaining results in reasonable agreement 
(in absolute units) with experimental data. They used 
the linearized time-dependent Ginzburg-Landau (TDGL) 
equation, which is identical to the Aslamazov-Larkin d(j| 
approximation for the microscopies. Linearized (Gaus- 
sian) TDGL is applicable far enough above the mean-field 
transition temperature Tj^^, where the order parameter 
fluctuations are small enough to neglect the nonlinear 
terms in the full Ginzburg-Landau free energy. To esti- 
mate the Nernst effect closer to Tc they also treated the 
nonlinearities in TDGL in the Hartree approximation . 

Prominent in the phase diagrams of the cuprate su- 
perconductors in a magnetic field is the vortex liquid 
regime, where short-range superconducting correlations 
are strong (because T < Tc), but thermal phase and vor- 
tex fluctuations disrupt the superconducting coherence 
on longer scales. The experimentallyobserved Nernst ef- 
fect remains large in this regime 0i M ^ IS j which is not 
accessible to analytic calculations due to its strong cor- 
relations and fluctuations. However, the TDGL equation 
with its full nonlinearity and driven by thermal fluctua- 
tions can be numerically simulated in this vortex liquid 
regime, and its Nernst effect measured and compared to 
experiments, as we demonstrate below. 

When the TDGL approximation is used to model prop- 
erties of cuprate superconductors, for which there is not a 
well-established more microscopic theory, or more gener- 
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ally when it is used well away from the zero-fleld critical 
temperature, as we do here, it constitutes a standard, 
widely-used, but essentially phenomcnological approxi- 
mation, ft has the virtue that most of the parameters 
that enter can be determined from experiments, as we 
discuss below. For BCS superconductors, the TDGL ap- 
proximation can be systematically derived from the mi- 
croscopies in some limited regimes [HI Il2j , but here we 
are instead using it more broadly, in regimes where its 
justiflcation is "only" phenomcnological. 

The Nernst effect in the vortex liquid can be described 
in terms of the vortices as the phase-slip voltage due to 
vortices being transported down the temperature gradi- 
ent as heat carriers. However, when one calculates the 
Nernst effect at higher temperatures in linearized TDGL 
, what is happening is that superconducting order pa- 
rameter fluctuations are transported along the temper- 
ature gradient and their phase patterns are twisted by 
this motion across the magnetic field, inducing a Nernst 
voltage; vortices play no role in the calculation. In a 
TDGL treatment of the vortex liquid regime, there is 
no clean distinction between order parameter fiuctua- 
tions and vortex motion, so both of these two different- 
sounding descriptions are in some sense correct. 

In this paper, we present the method for and results 
of simulations of thermoelectric transport in type-II su- 
perconductors in the vortex liquid regime. We simulate 
the TDGL equation with thermal noise. We work in the 
strongly type-II limit, n ^ 1, where the magnetic field in 
the sample is assumed to be uniform and not fluctuating. 
Mostly we work in two dimensions, but we also examine 
the crossover from two-dimensional to three-dimensional 
behavior. In our study of interlayer couplings we see 
some indication that a substantial part of the entropy 
carried by the vortices in the vortex liquid may be the 
configurational entropy of their positions. 

Our simulations involve a dimensionless tunable pa- 
rameter (77) which sets the strength of the thermal fluc- 
tuations in the sample and hence can range from nearly 
mean-field behavior with very weak fiuctuations to very 
strong fluctuations with strongly suppressed supercon- 
ductivity. We show that for an intermediate value of this 
parameter our results are in reasonably good qualitative 
and quantitative agreement with available experimental 
Nernst-effect data on overdoped La2_xSra;Cu04 (LSCO) 
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In this comparison all the other parameters in the 
simulations are set by experimentally measured quanti- 
ties. We also discuss in brief the results of our simulations 
to model experiments on Bi2Sr2CaCu208+rr (BSCCO) 



II. SOME THERMO-ELECTRO-MAGNETIC 
TRANSPORT THEORY 

We briefly review some of the general theory of thermal 
and electrical transport coefficients in linear response in 
the presence of a magnetic field. Consider a sample at 
temperature T subjected to small gradients in the poten- 
tial V(/) and temperature VT. One can then very gen- 
erally write down the transport current densities in the 
system to linear order in V(/) and VT as 



aV0 - kVT 



(1) 



(2) 



where J^j. is the charge transport current density and J^. 
is the heat transport current density, a, a, a and k are 
the electrical conductivity, thermoelectric, electrother- 
mal and thermal conductivity tensors. The Onsager re- 
lations for the transport coefficients tell us that a = Ta. 
The Nernst coefficent {v) can be defined in a configura- 
tion with a magnetic field B|lz and VT||x along with the 



condition Jf 



as 
E,, 



1 ^xy^xx ^xx^x 



BVT B 



(3) 



'xy 



If the system shows no Hall effect (like the one we have 
simulated), then a^y = and 



Ba,, 



(4) 



One way to obtain transport coefficients from a simu- 
lation is via the Kubo relations, thus measuring current- 
current correlations in the equilibrium dynamics. Alter- 
natively, one can apply an electric field and/or a temper- 
ature gradient and measure the resulting currents. We 
have used both methods, finding that our implementation 
of the latter method gives higher signal-to-noise ratio per 
unit computer time in the regime we have been studying. 
The currents flowing in the sample have magnetization 
parts in addition to the transport parts considered above. 
One could write down linear response relations for the to- 
tal current densities analogous to the ones for the trans- 
port current densities but the coefficients appearing in 
those would not obey the Onsager relations. The total 
currents in the system can in general be computed in a 
straightforward manner but one has to extract the trans- 
port parts from these to calculate the transport coeffi- 
cients. This involves identifying and subtracting out the 



magnetization currents. The procedure we use to that 
end in this work is outlined below and is based on the 
arguments of Cooper et al. 13] . 

Let Jtot('") ^iid Jtot(i") be the total charge and heat 
current densities (transport -I- magnetization) at any 
point r in the sample. If 0(r) is the electric potential 
at r, there exists a total energy current J^^. (r) such that 



J,«„,(r)=jf„,(r)-0(r)JLt(r) 



(5) 



A similar relation holds between the transport parts of 
these current densities. 



jg(r) = J,^(r)-0(r)J^,(r) 



(6) 



The average electric current density Jj?j, and heat current 
density are given by averaging J^^ (r) and (r) over 
the whole sample. For a two dimensional sample 



and 



3t, = ^ J^JtAr)dS 



(7) 



(8) 



where S is the area of the sample. Similar relations also 
hold for three dimensional samples. From Eqns. (6) and 
(8), we obtain 

= ^ (/^ 3fAr)dS - ^ 4>{r)3Ur)dS^ . (9) 



We know that 



JtAr) = J^,(r)+J 
3g,{r) = Jt^(r)+J 



ma; 
E 

mag 



(r) 
(r) 



(10) 



where JJjjag(r) and Jmag(i") ^i'^ the charge and energy 
magnetization current densities. It can be shown on gen- 
eral grounds that there exist charge and energy mag- 
netization densities M'^(r) and M^(r) such that 



mag 
E 

ma; 



(r) = V X M'=(r) 
(r) = V X M^(r) 



(11) 



The charge magnetization of the system is nothing but 
the conventional magnetic moment density M(r). If the 
material surrounding the sample is assumed to be non- 
magnetic, both M'^(r) and M^(r) vanish outside the 
sample. The magnetization currents are mostly at the 
boundaries because that is where the change in magne- 
tization is the largest. Thus we see that 



3Ur)dS 
3^,.ir)dS 



jfot(r)d^ 



(12) 
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From Eqns (9) and (12), we obtain 

3?, = ^ (^J^3Ur)dS - J^ cj,{r)JUr)dS^ . (13) 

It should be noted that there is no "heat magnetiza- 
tion density" M'3(r) analogous to M''{r) and M^(r) 
for which the heat magnetization current JmagC'") = 
V X M'3(r). A consequence of this is that unlike and 



Is JS 
In fact, it can be shown that 



Js 

shown that 
jgaK(r) = V0 X M(r) . 



(14) 



III. SIMULATIONS 



We measure the transport coefficient ctxy by turning 
on an electric field — Vf/fHy and a magnetic field B||z and 

-Q(x) 

measuring J^^ , the x component of the average heat 
transport current density at constant temperature. This 
is done in a system with mixed boundary conditions: it 
has free surfaces at the boundaries normal to y, while it 
has periodic boundary conditions along the x-direction 
(and along z when we study 3D). If, as in our model, the 
system does not have a Hall response, under these con- 
ditions the X component of the electric transport current 

vanishes at all points, J^j'^\r) = 0. Thus from Eqns. (2) 
and (13), we obtain 



(X 



(15) 



Similarly, Uxy can be measured by introducing a tem- 
perature gradient VT||y and magnetic field B||z with no 

external electric potential and measuring J[j. 
Eqns. (1) and (12), 



■e(x) 



From 
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r)dS 



VT 



(16) 



The two coefficients axy are axy can thus be obtained by 
measuring only the total currents in the system. They 
are related to each other by the Onsager relation and to 
the Nernst coefficient by Eqn. (4). 

The TDGL equation that we have simulated is 



4' 



(17) 



2m 

-ao{T - T*^^)^ - b\^\H + C(r, t) 



We take r to be real so there is neither Hall effect nor 
Seebeck effect f6| . This approximation of leaving out the 



Hall effect should be reasonable as long as the Hall angle 
is small, as it generally is in the vortex liquid regime 
of the cuprate superconductors. We work in the type-II 
limit where the magnetic field is assumed to be uniform 
and not fiuctuating. The noise correlator is 

< C (!■', t')C(r, t) >= 2TkBT6{v - r')S{t - t') . (18) 

The current operators we need are 0] 



-I (** 

2m* 

~2m* ^~dt 



V - 
V 



h 

ie' 
H 



A)*)+c.c. (19) 
A I + c.c. 



The TDGL equation is the simplest dynamical stochas- 
tic equation one can obtain from the Ginzburg-Landau 
free energy functional. As a consequence of its simple 
relaxational dynamics, it does not explicitly conserve ei- 
ther total charge or total energy. It is implicitly in con- 
tact with and exchanging energy and charge with lo- 
cal reservoirs. In reality these reservoirs are presumably 
the phonons, quasiparticles and other excitations not in- 
cluded in the TDGL equation. Nevertheless, the charge 
and energy currents carried by the superconducting or- 
der parameter 5* and its fluctuations are as given above 
and can be measured and their transport properties de- 
termined within this model. Charge and energy conser- 
vation can explicitly be taken into account in microscopic 
derivations of the TDGL theory, where only the contri- 
bution due to superconducting fluctuations is retained. 
It can however be shown that other effects (normal state 
contributions etc.) do not contribute as significantly to 
the Nernst effect as superconducting fluctuations 

An important feature of TDGL for the present work is 
that the parameter r, which sets the time scale, does not 
enter in the values of the transverse thermoelectric coef- 
ficients that we are studying, axy and oixy This allows 
quantitative comparison to experiments to be done with- 
out estimating r, which is fortunate, because the value 
of r is quite uncertain. 

This TDGL equation with noise has a dimensionless 
parameter that gives the strength of the thermal fiuc- 
tuations. To remove all the dimensional quantities, we 
use the zero temperature coherence length, as our 
unit of length; T*^^, fcsTj^^, A and r/iaoTf^) as our 
units of temperature, energy, magnetic flux and time, 
respectively; and the zero-temperature order param- 
eter magnitude, as our unit for the order parameter. The 
resulting TDGL equation is then 

(dt -f i^) * ( V - iAf - (T - 1)^- - l^p^f + C(r, t) , 

(20) 



with noise correlator 



< C(r', t'Xir, t) >= 2T]TS{r - r')S{t - t') 



and 



MF 



(21) 



(22) 
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is the fluctuation parameter, where d is the number of 
spatial dimensions. When 77 ^ 1, the actual critical tem- 
perature Tc is near the mean-field critical temperature 
T^'^, while if 77 > 1 then < T^^. The model has 
four parameters that remain after scaling to the units 
specified above: the temperature, the magnetic field, the 
cutoff, and the strength 77 of the noise. 

We initially simulate a two-dimensional system, thus 
ignoring interlayer coupling. We discretize space and 
time. The spatial grid spacing is taken to be ^o- This 
rather coarse spatial grid is used to minimize the com- 
puter time needed to simulate samples large compared to 
this microscopic length. There are noticeable quantita- 
tive effects of using such a coarse grid. For example, in 
our units Hc2{T = 0) = 1.18 with this grid, while it is 1.0 
for the continuum model. However, recognizing that the 
model we are simulating is rather simple, so should not 
be expected to give highly precise quantitative results, 
our goal in the present study is to be quantitatively only 
as accurate as might be expected of such a simple model. 
Although we do not know how accurate that really is, 
the general precision to which we have chosen to work 
is roughly 10% to 20%. The quantitative shifts due to 
our using a coarse grid in the simulations are at most of 
this magnitude. We have also verified that using a finer 
grid spacing would not alter our results in any significant 
way. Reduction of the grid spacing adds more degrees 
of freedom, resulting in increased fluctuations, and sup- 
pression of Tc and Hc2- However, if one also lowers the 
value of the fluctuation parameter 77 a little, this increase 
in the fluctuations can be removed, and the resulting be- 
havior of the system is only weakly dependent on the grid 
spacing. 

In the simulations, the time step is chosen to be be- 
tween 0.02 and 0.1 in our units, with shorter steps at 
higher temperatures; we check that our results are not 
affected by doubling the time step. A first-order Eu- 
ler method is used, with 10000 to 20000 time steps 
used for equilibration, and 30000 to 50000 steps for av- 
eraging. The two-dimensional data we show are for a 
50x50 square grid. We use the gauge- invariance of the 
TDGL equation to work as much as possible in terms 
of gauge- invariant quantities like \^\ and the gauge- 
invariant phase differences 

uj{r,r',t) ^9{r',t) -9{r,t) - A • dl , (23) 

where 4* = |4'|e**. The current densities are also written 
in terms of these quantities. 

At _ff = 0, for the two-dimensional system there is a 
Kosterlitz-Thouless transition at Tc, which we locate 
by measuring the helicity modulus. 

To obtain dxy, we measure the transverse energy cur- 
rent due an electric field. We then use the Onsager re- 
lation to obtain axy — axy/T. We find that the sig- 
nal to noise ratio for this measurement is significantly 
smaller than in a direct measurement of a^y We have 
numerically verified the validity of this Onsager relation 



by measuring a^y both ways. This also serves as a check 
that we are indeed measuring the proper transport co- 
efficients. We also confirm that our results agree with 
the analytic results 7| in the higher-temperature regime 
where linearized TDGL applies. 

The temperature enters the TDGL equation both in 
the term linear in the order parameter and in the inten- 
sity of the noise. It can be shown that introducing a gra- 
dient in the linear term only adds to the magnetization 
currents, and thus for measuring the transport currents 
it is sufficient to include the temperature gradient only 
in the noise term. We have also verified this numerically. 



IV. RESULTS 

Some results for a^y from the two-dimensional simu- 
lations are shown in Fig. 1. We have chosen to compare 
these results to the experimental results of Wang et al. 
1^ obtained from Nernst and resistivity measurements on 
an overdoped LSCO sample with x = 0.2 and Tc = 28K. 
For this comparison, we present the simulation results 
in SI units using the parameters from the experiment, 
Hc2{T = 0) 45T (thus ^ ^ 27 A) and layer spacing 
s = 6.6A. The fluctuation parameter has been adjusted 
to 7y = 0.42 to give reasonable quantitative agreement 
between simulation and experiment. We set the temper- 
ature scale using T*^^ = 40K, which gives Tc = 28K in 
the simulation. Note that axy is not directly measured 
in the experiments. For a system with negligible Hall 
effect, axy can be obtained from the measured Nernst 
effect and the longitudinal resistivity. Thus for this com- 
parison we require both these measurements to be made 
on the same sample over a substantial portion of the vor- 
tex liquid regime. This requirement seriously limits the 
number of experimental results we can compare to. 

Looking at Fig. 1, we see that there is reasonably good 
general qualitative and quantitative correspondence be- 
tween the simulation and experiment, considering how 
simple the model is and that we have adjusted only the 
strength of the thermal fluctuations 7/ to get this agree- 
ment. The other parameters are all dictated by experi- 
ment. Our results are thus consistent with the proposi- 
tion that the superconducting fluctuations modelled by 
the TDGL equation produce most, if not all, of the con- 
tributions to the large Nernst effect seen in the vortex 
liquid regime for this cuprate sample. Of course, the 
agreement between experiment and this simple model is 
not perfect. One difference is seen at high field in the 
higher-T panel of Fig. 1, where the simulations give a 
nearly-constant axy, while the experiment shows more 
variation. 

One effect that is not included in our two-dimensional 
simulations is interlayer coupling. For the cuprates, this 
can become important in the vortex liquid at low fields 
in the vicinity of Tc, where there is a crossover to three- 
dimensional behavior as the correlation length normal 
to the layers becomes larger than one layer spacing. 
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Note that the difference between simulation and exper- 
iment is large at low fields and T = 24K. We consid- 
ered the possibility that this difference might be due to 
the crossover to three-dimensional behavior, which is not 
present in our two-dimensional simulation. To investi- 
gate this crossover, we have also simulated the Lawrence- 
Doniach version of the TDGL equation. In our units 
this is 



20 
H(T) 




FIG. 1: Experimental data obtained from measurements of 
the Nernst effect on overdoped LSCO {x — 0.2) along with the 
results of our two-dimensional simulations. The parameters 
used in the simulation are: = 40K, Hc2{0) = 45T, s = 

6.6A, rj = 0.42. The critical temperature is Tc ~ 28K in both 
simulation and experiment here. 



J 



(dt + = (Vi - tAi_)Hj -{T- 1)*, + J(e 



-isA, 



J + 1 



(24) 



r 



where J is the interlayer Josephson coupling per unit 
area, and j is the layer index. The transition temperature 
is located by finding the intersections of the fourth-order 
cumulant (Binder ratio) curves (|'0|'')/(|V'P)^for different 
system sizes as functions of temperature at zero magnetic 
field. 



Fig. 2 shows results for function of field 

for various values of interlayer coupling J at tempera- 
ture T — 24K, from samples of size 15x15x10 layers. 
For each J we set the the temperature scale {T^^) so 
that Tc = 28K as in the experiment. The values of 



the other parameters in the simulations are as used in 
Fig. 1. As expected, the the largest value of J shows 
the largest deviation from the two-dimensional behavior. 
Somewhat better correspondence between the behavior 
of ot'xy obtained in simulation and experiment is obtained 
by adding this interlayer coupling. 

An interesting point to note in Fig. 2 is that at low 
fields a rather small interlayer coupling J produces a 
strong decrease in axy Describing this in terms of vor- 
tices, as seems appropriate in this low-field regime below 
Tc, the force per unit length on a vortex line due to the 
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FIG. 3: (Left) Data from the simulation with s = 15.35A, 
Tc = 90 K, Hc2 = 160 T and 77 = 0.475. (Right) Experimen- 
tally measured 0] value of obtained from Nernst measure- 
ments on optimally doped BSSCO (2212) with Tc= 90 K and 
Hc2{0) ^ 160 T. 
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H(T) 

FIG. 2: Three-dimensional simulation results for T — 24K, 
with interlayer coupling J. Except for T^'^ , the parameters 
used here are the same as in Fig. 1. See text for more details. 

temperature gradient is proportional to axy- Naively, 
this force is also set by the entropy per unit length trans- 
ported by a moving vortex. This transport entropy may 
come from internal degrees of freedom within one vor- 
tex, or it may come from the configurational entropy of 



the many possible spatial arrangements of the vortices. 
The latter does not involve internal excitations within 
the vortices. In the quasi-two-dimensional layered super- 
conductors that we are modelling, a vortex line running 
normal to the layers consists of "pancake vortices" in each 
layer. A weak interlayer coupling has little effect on the 
internal properties of a single pancake vortex, but it does 
produce an attraction between vortices in adjacent layers 
that reduces their relative motion, and thus can substan- 
tially reduce the configurational entropy of the vortices. 
The fact that at low fields a weak interlayer coupling 
greatly reduces the transport entropy of the vortices, as 
indicated by the reduction of Uxy seen in Fig. 2, thus sug- 
gests that in this regime the configurational entropy of 
the vortices is a large part of their transport entropy, at 
least within this layered TDGL model. 

We have also performed simulations to model the 
BSCCO sample studied by Ri et al. 0]. The experi- 
mental data and the results of the simulation are shown 
in Fig. 3. The quantity plotted is the measured "trans- 
port energy" of the vortices, defined as = TcjiQUxy 
We used a two-dimensional simulation with Tc ^ 90 K, 
Hc2 = 160 T, s = 15.35A and rj = 0.475. Again, there 
is reasonable quantitative agreement, in absolute units, 
between the experimental data and the simulations. The 
clearest deviation here occurs for T < 60K and H > 6T, 



where the sign of the temperature derivative of Urj, ap- 



pears to differ between simulation and experiment. 
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In conclusion, we have measured the transverse ther- 
moelectric transport coefficient a^y in simulations of the 
vortex-liquid regime of superconductors modeled by the 
TDGL equation with thermal fluctuations. Wc find that 
our simulations of a two-dimensional superconductor re- 
produce reasonably well much of the qualitative and 
quantitative features of available experimental data from 
some cuprate high-temperature superconductors. We 
have also studied the crossover from two-dimensional to 
three-dimensional behavior in a layered superconductor. 
This crossover study indicates that the configurational 
entropy of the vortices may constitute a large part of 



their transport entropy at low magnetic field just below 

c- 
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